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ABSTRACT 

We analyse the skewness of the hne-of-sight velocity distributions in model ellipti- 
cal galaxies built through collisionless galaxy mergers. We build the models using large 
N-body simulations of mergers between either two spiral or two elliptical galaxies. Our 
aim is to investigate whether the observed ranges of skewness coefficient (/13) and the 
rotational support (y/cr), as well as the anticorrelation between and V, may be 
reproduced through collisionless mergers. Previous attempts using N-body simulations 
failed to reach V/a « 1 — 2 and corresponding high /13 values, which suggested that gas 
dynamics and ensuing star formation might be needed in order to explain the skewness 
properties of ellipticals through mergers. Here we show that high V/a and high h^, are 
reproduced in collisionless spiral-spiral mergers whenever a central bulge allows the 
discs to retain some of their original angular momentum during the merger. We also 
show that elliptical-elliptical mergers, unless merging from a high-angular momentum 
orbit, reproduce the strong skewness observed in non-rotating, giant, boxy ellipticals. 
The behaviour of the /13 coefficient therefore associates rapidly-rotating disky ellip- 
ticals to disc-disc mergers, and associates boxy, slowly-rotating giant ellipticals to 
elliptical-elliptical mergers, a framework generally consistent with the expectations of 
hierarchical galaxy formation. 
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numerical simulation 



1 INTRODUCTION 

Elliptical galaxies are frequently divided into two sets with 
different characteristics. High luminosity ellipticals showing 
boxy isophotes are mainly supported by velocity dispersion, 
and they are commonly associated to extended X-ray haloes 
and are radio loud. Low luminosity elliptical galaxies present 
a high amount of rotation being rotationally supported. 
They show disky deviations in their elliptical isophotes and 
usually they are not associated to X-ray emission nor to ra- 
dio loud objects (Lauer 1985, Bender 1988, Bender et al. 
1989). 

Explaining such dichotomy is a basic requirement of el- 
liptical galaxy formation theories. Early N-body work, which 
pointed out that 1:1 collisionless disc galaxy mergers give as 
a result a body that resembles a boxy, slow-rotator elliptical 
galaxy, while mergers with mass ratios of the order 3:1 give 
disky ellipticals (Barnes 1998; Naab et al. 1999; Bendo & 
Barnes 2000, hereafter BBOO; Naab & Burkert 2001, 2003) 
ignored the basic fact that the two classes of ellipticals are 
mainly segregated by luminosity, and hence, mass. In hi- 
erarchical scenarios for galaxy growth, because successive 



merging moves galaxies into progressively bulge-dominated 
systems, elliptical galaxies are end-points, hence they domi- 
nate the high end of the mass distribution. In these schemes, 
massive z = Q ellipticals, would be expected to form through 
mergers of ellipticals. Because of this, explanations for the 
dichotomy in elliptical galaxy properties in a hierarchical 
context will naturally be related to differences between the 
galaxy types of the precursors. 

In this letter, we elaborate on the types of precur- 
sors that may have merged to form elliptical galaxies by 
analysing the line-of-sight velocity distributions (LOSVD) 
in simulated N-body merger remnants. LOSVDs of ellip- 
tical galaxies may be obtained from high signal-to-noise 
(S/N) spectroscopy, and are customarily modeled by fitting 
a Gauss-Hermite series (e.g., van der Marel & Franx 1993). 
Two indices, hs and hi, parametrize the strength of asym- 
metric and symmetric deviations from Gaussian shape, re- 
spectively. Bender, Saglia & Gerhard (1994, hereafter BSG) 
show that /i3, the skewness of the LOSVD, has opposite sign 
to the rotation velocity: the LOSVD is steep in the leading 
side of the rotation, and has a more extended tail on the 
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trailing side. Furthermore, BSG show that for eUipticals as 
a class, fta and the rotational support parameter V/a follow 
a distinct pattern. We show their main result (data kindly 
provided by R. Saglia) in Figure^ which plots the skewness 
coefficient /13 against V/a, for a number of galaxies; each 
point corresponds to a local position along a major-axis slit, 
and all galaxies are shown in the same diagram. Clearly, 
/i3 and V/a have opposite signs. At low |V/o-|, which corre- 
sponds to giant ellipticals and to the central parts of other 
ellipticals, \hz\ increases rapidly with |V^/o-|. At high |V^/o-|, 
which maps intermediate and outer radii in intermediate- 
luminosity ellipticals, /13 reaches up to |/i3| ~ 0.2. 

Several works have studied what types of LOSVD are 
generated by mergers. N-body studies are useful tools as 
they allow to dissect the remnant LOSVDs into components 
by, e.g., the contribution by each merger precursor, or by the 
contribution of each orbit type. The first study of LOSVDs 
of N-body merger remnants (Balcells 1991) showed that 
mergers of unequal ellipticals lead to asymmetric LOSVDs 
similar to those observed in kinematically decoupled cores 
(steep leading sides) and that the asymmetry is not due to 
an accreted disc but to dynamical transformations in the 
primary as a result of the merger. Bendo & Barnes (2000, 
hereafter BBOO) pioneered the use of Gauss-Hermite analysis 
of LOSVDs in N-body merger remnants, and demonstrated 
that the shape of the LOSVD is intimately linked to the 
orbital makeup of the remnant. In short, box orbits have 
no net rotation, hence their contribution to the LOSVD is 
fairly symmetric and centered on V = Q. Short-axis tubes 
show net rotation in merger remnants, hence their contri- 
bution to the LOSVD peaks at F > 0. Whenever short-axis 
tube orbits dominate over box orbits, the global LOSVD 
will be steep on the leading side (leading to LOSVD similar 
to those observed by BSG), while dominance of box-orbits 
leads to a peak at V ~ and an extended tail on the leading 
side. 

Naab & Bukert (2001, hereafter NBOl) were the first 
to point out that collisionless major mergers of disc galax- 
ies (mass ratios 1:1 and 3:1; disc: bulge ratio 3:1) lead to 
LOSVDs with opposite behaviour to observations: has 
the same sign as V/o. They noted that their models can 
be made to agree with the observational trends (opposite 
signs for and V/a) by postulating that ellipticals harbor 
large-scale stellar discs (15% in mass, scale-length h ~ Rcs) 
formed after the merger from left-over gas. The statistics of 
isophotal diskiness suggest that such discs might in fact ex- 
ist in most intermediate-luminosity ellipticals (Rix & White 
1990). 

The dissection of LOSVDs into their orbit constituents 
discussed above clarifies why collisionless major mergers lead 
to the wrong correlation between hz and V/o: such merger 
models have too many box orbits. The discs postulated by 
NBOl simply add short-axis tubes and compensate for the 
excess of box orbits. It is useful to enquire about ways to 
increase the relative proportion of z-tube orbits which are 
more directly related to the merger dynamics. Box orbits 
are favored by cores with shallow density profiles, hence the 
limited spatial resolution of the simulations might lead to 
too many box orbits. Limited resolution of the central parts 
also causes an earlier disruption of the merging galaxies, 
and to a remnant which is less supported by net rotation, 
again enhancing box orbits. Recent results on collisionless 
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Figure 1. /13 vs. V/a for 44 near-by elliptical galaxies from BSG. 



mergers of discs (Gonzalez-Garci'a & Balcells 2005, hereafter 
GGB05) show that the bulge component helps stabilizing 
the disc during the merger. Right before merging, the pro- 
genitor discs are still largely in place, delivering more an- 
gular momentum to the remnant. That could enhance the 
proportion of short-axis tube orbits in the remmnant. Fi- 
nally, it should be noted that, in the study of this question, 
mergers between elliptical galaxies have not been considered 
so far. 

In the present letter we consider merger simulations in- 
volving two disc galaxies with large bulges or two elliptical 
galaxies, varying the initial orbital parameters of the en- 
counters. In Section |5| we present our simulations, in Sec- 
tion 121 we explain the LOSVD derivation procedure. Sec- 
tion |1| presents our results and finally in Section |3 we dis- 
cuss the implications from our results for the current view 
of elliptical galaxy formation. 



2 MODELS 

The disc galaxy models comprise a halo, a disc and a bulge 
built following Kuijken & Dubinski (1995). Masses are in 
the ratios halo:disc:bulge = 36:5:1 and 19:2:1; other model 
parameters are those of the dbh models of GGB05, but with 
a larger number of particles. We note that the initial mass 
ratio between disc and bulge is 5:1 for one of the models 
and 2:1 for the other, this is a somewhat larger bulge than 
used by Naab & Burkert (2001, 2003). Models with different 
masses are scaled up version of low mass systems following 
a TuUy-Fisher like relation. 

The elliptical galaxy models are Jaffe (1983) non- 
rotating isotropic models without a dark matter halo. For 
details, see Gonzalez- Garcia & van Albada 2005a. The Jaffe 
scale parameter was set to rj = 1; truncation at r = 10 x rj 
means that the true half-mass radius is ri/2 = 0.82. Systems 
of masses different from one are scaled up versions of M = 1 
systems following a Fish (1964) law. The initial model was 
let to evolve in isolation for 10 crossing times. After a mild 
re-organization, the system reaches equilibrium. This equi- 
librium configuration is the one used in the merger models. 

Model units are Mdisc = 1, /idisc ~ 1 and G = 1. A set 
of physical units that match the dbh models to the Milky 
Way are: [M] = 3.24 x 10^^ Mq; [L] = 14.0 kpc; [T] = 
4.71 X 10^ yr, with [v] = 315 km/s. For scaling our Jaffe 
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Table 1. Initial configurations for the models. (1) Model name, 

(2) type of merger, disc-discs (S+S) or elliptical-elliptical (E-l-E), 

(3) Mass ratio, (4) initial orientation of the spin, (5) initial sepa- 
ration, (6) orbit ellipticity and (7) pericenter distance. The latter 
two parameters refer to the Keplerian orbit of same initial condi- 
tions. 
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models to elliptical galaxies the set of units would be: [M] — 
Mj = 4 X 10^^ Mq; [L] = rj = 10 kpc; [T] = 2.4 x 10^ yr, 
and finally [v] = 414 km/s. 

Parameters for the mergers analysed here are given in 
Table All experiments involve direct orbits. We include 
four spiral-spiral mergers, and four elliptical-elliptical merg- 
ers. For the discs, we chose slightly sub-parabolic orbits to 
accelerate merging. Remnants of S+S mergers with mass ra- 
tios 4:1 to 7:1 are reported to have characteristics between 
those of ellipticals and spirals (Naab & Burkert 2003, Bour- 
naud et al. 2004) thus a choice of mass ratios 1:1 to 3:1 
seems adequate for this study. For the ellipticals, merger 
orbits were parabolic. Two equal-mass E+E mergers with 
different orbital angular momentum were computed, to in- 
vestigate the dependence of the LOSVDs on the orbital an- 
gular momentum. 

For the disc models we use 2.75 x 10"" particles for the 
halo, 1.75 X 10^ particles for the disc and 9 x 10* for the 
bulge. This makes a total of 10^ particles for the merger rem- 
nant. Each elliptical galaxy is modeled with 10^ particles. 
We run the simulations with the parallel version of GAD- 
GET (Springel et al. 2001, Springel 2005) on the TAG Be- 
owulf Gluster with 16 GPU's and on the MareNostrum clus- 
ter of the Barcelona Supercomputer Genter^. Energy con- 
servation is good, and variations are kept below the 0.5%. 
Independent softening is used for each particle type. The 
softening is taken as one fifth of the half mass radius of the 
component. After merger is completed we let the remnants 
to evolve in isolation for 10 half-mass radius crossing times 
to ensure that the inner parts are close to virial equilibrium. 



3 LOSVD DERIVATION 

We obtain LOSVDs for the merger remnants in the following 
way. We choose a point of view at random and project the 
particle distribution. We then derive a surface density plot 
by binning the space and counting the number of projected 
particles in each bin. Finally we define isodensity contours 
and fit ellipses deriving values for the ellipticity, position 

^ http://www.bsc.es 



Figure 2. h-j vs. V/a for the disc mergers. Top row models with 
equal mass mergers, bottom row 3:1 merger. Left panels bulge- 
to-disc ratio 1:5, right models with B:D of 1:2 (a) Model Sdll. 
(b) Model 2dll. (c) Model 5d31. (d) Model 2d31. 

angle and the 04 parameter. Then we place a slit along the 
major axis of our ellipses up to one effective radius i?cff- 
We have binned the slit in ten spatial bins and the velocity 
interval in 50 bins. We find the radial projected velocity 
and the number of particles in each bin in velocity for each 
bin in the slit. In this way we obtain a line-of-sight velocity 
distribution. 

Finally we fit the LOSVD by a Gaussian and the resid- 
uals by a Gauss-Hermite polynomial as given by van der 
Marel & Franx (1993) and BSG. We repeat this process for 
each remnant for 90 randomly chosen points of view. 

From the fitting procedure we obtain a value for V , the 
velocity centroid of the distribution at each spatial bin along 
the slit, the velocity dispersion a, and the amplitude of the 
third Hermite polynomial ^3, which is a measurment of the 
skewness of the distribution. 



4 RESULTS 

Following BSG, we plot together /i3 vs. V/a for all spa- 
tial bins of all slits for the merger remnants. LOSVDs are 
mapped out to 1 Rca- We first analyze the remnants of 
disc galaxy mergers, see Figure |5| Galaxy mass ratios, and 
bulge-to-disc ratio (B:D or D/B in the following) for the 
precursors, are listed in the legends. For equal-mass mergers 
(Fig. 121 top panels), /13 vs. V/a markedly depends on the 
precursors' D/B. A small bulge (D/B=5, Fig. |2K) leads to 
low \V/(t\ < 0.5, and small asymmetries, \hz\ < 0.03, with 
hz positively correlated with V/o. Such positive correlation, 
opposite to that observed in real ellipticals, was also found 
by NBOl, and would easily lead one to conclude that coUi- 
sionless, equal-mass mergers of spirals yield wrong LOSVDs. 
As it turns out, the LOSVDs depend on other model param- 
eters besides the galaxy mass ratio. Figure |2J3 corresponds 
to the same equal-mass merger as in Figure |2^ but with a 
larger bulge (D/B=2). This model shows higher rotational 
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Figure 4. Combined /13 vs, y/cr relation for S+S and E+E mod- 
els. 



Figure 3. /13 vs. y/cr relation for the four E+E merger models. 

(a) Model ell with mass ratio 1:1 and small impact parameter, 

(b) Model e31 with mass ratio 3:1, (c) Model e51 with mass ratio 
5:1 and (d) Model elld with mass ratio 1:1 and a large impact 
parameter 



support within 1 RcS, jV'/c"! < 1-3, and higher asymmetries, 
l/isl < 0,1, More important, is anticorrelated with V/o, 
except for a small region at \V/(t\ ~ 0, We conclude, in agree- 
ment with NBOl, that equal-mass collisionless mergers of 
spirals do not reproduce the observed LOSVD asymmetries 
of elliptical galaxies. But note that a significant amount of 
hz-V/o anticorrelation can be obtained through these merg- 
ers, if the precursor host massive bulges. 

Figure shows the results for our 3:1 merger remnant 
with D/B=5, We find larger values for V/cj{e (-2,2)); it 
is well known that 3:1 S-j-S mergers lead to remnants with 
higher rotation (BBOO; GGB05), The /13 asymmetry term 
has opposite sign to V/o, for high V/a, whereas it has the 
same sign as V/a for low V/a. The /13 - V/a distribution 
in Figure |2; is quite unlike the observations. As was the 
case with the equal-mass mergers, this behaviour is sensi- 
tive to D/B, Figure |5Ji corresponds to the same model as in 
Figure except that now D/B=2, Here the and V/a 
are anticorrelated for all values of V/a. While the h-j, -V/a 
distribution of Figure|5Ji fails to match entirely the BSG ob- 
servations, we highlight that our model reproduces the cor- 
rect sign for the hz -V/a correlation. Therefore, it is is not 
true that collisionless mergers lead to positively-correlated 
/13 and V/a. 

We argued in §Qthat the hz-V/a distribution is mostly 
given by the relative weights of box and short-axis tube or- 
bits. Matching observations (i,e,, obtaining anti-correlated 
hz and V/a) requires increasing the fraction of short-axis 
tubes. It is likely that our D/B=2 models reach high hz 
and high V/a because they have deeper central potentials 
than those of NBOl, owing to a more massive bulge. Deeper 
galaxy central potentials stabilize the merging discs against 
bar distortions during orbital decay. Spin angular momen- 
tum is then deposited deeper into the remnant's luminous 
body, which should naturally lead to more z-tube orbits. In 
contrast, discs which bar-distort during orbital decay trans- 



fer their spin angular momentum to the haloes, leading to 
tumbling prolate or triaxial remants with little net rotation 
(GGB05), Differences in the cuspiness of the halo potentials 
may also lead to different relative proportions of box and 
z-tube orbits in ours and NBOl models. 

Two conclusions arise from the previous analysis. First, 
that the hz-V/a distributions from S-j-S mergers depend on 
subtleties of the models such as the D/B or the shape of 
the inner halo potential; and, second, both ours and NBOl 
models suggest that the full hz-V/a distribution observed in 
elliptical galaxies (our Fig, cannot be reproduced by S-j-S 
mergers alone. 

Figure |3 shows the results for our mergers of ellipti- 
cal galaxies. The low-impact parameter mergers (Fig, |3Ji.bc) 
lead to similar distributions irrespective of the galaxy mass 
ratio: low V/a, large range of hz, and hz-V/a anticorrela- 
tion. These models occupy the same region of the hz-V/a 
plane as the giant, boxy, slowly-rotating ellipticals of BSG, 
Model elld (Fig,|^) slightly deviates from this behaviour; 
V/a reaches higher values, and hz is overall smaller, and 
slightly uncorrelated with V/a. 

Thus, a way to populate the region of low V/a that 
should not be overlooked is through mergers of E's, The 
outcome of these mergers show a range of properties, in- 
cluding diskiness (Gonzalez- Garcia & van Albada 2005a, 
2005b); model elld shows that E-fE mergers may even lead 
to large values of V/a although we must note that out ellip- 
tical galaxies models do not have a dark matter halo which 
would lead to a smaller transfer of angular mometum to the 
luminous parts (Gonzalez-Garcia & van Albada 2005b), 

What would the hz-V/a look like for a population of el- 
liptical galaxies formed through mergers with varying mass 
ratios and types of the precursors? In Figure 2] we simply 
overplot the hz-V/a distributions from all the models. The 
combined distribution captures some of the most salient 
features of the BSG hz-V/a distribution of real galaxies 
(Fig, 0, namely, a cloud at low V/a with high hz, and 
two wings at high- V/a, with anticorrelated hz- Only one 
discrepancy remains, our models yield many points of view 
with intermediate V/a and hz ~ 0, which are not observed 
in the BSG figure. 
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5 DISCUSSION 

Figure |l]suggests that the BSG h^-V/a distribution may be 
explained by S+S mergers leading to the rapidly-rotating 
ellipticals, and by E+E mergers leading to giant, slowly- 
rotating ellipticals. There is no single coUisionless origin to 
explain the entire BSG correlation. But coUisionless mergers, 
when including E-l-E mergers as well as S-l-S mergers with 
bulges as massive as in our models (B/D=l:2, i.e., similar 
to the Galaxy), explain the main features of the hz-V/o 
distribution, providing an interesting additional mechanism 
to that based on the NBOl assumption that a stellar disc 
must have formed after the merger from left-over gas. 

The h^-V/a relation may be seen as the imprint in the 
LOSVD of the angular momentum present, either of spin or 
of orbital origin, in the systems that merged to form present 
day ellipticals. Hence, LOSVDs give us information about 
the types of the progenitors and the orbits they were placed 
on. 

An interesting corollary of our explanation involving 
S/E precursors for the high/low V/a domains is that giant, 
slowly-rotating, boxy ellipticals, which dominantly reside at 
low V/a, must have formed through mergers of ellipticals. 
While, at low redshifts, intermediate-mass (4:1 - 10:1) merg- 
ers are more likely than major E-l-E mergers (Bournaud 
et al. 2005), observational evidence for the importance of 
E-l-E mergers at higher cosmological distances is mounting 
(van Dokkum et al. 2003, BeU et al. 2004, Tran et al. 2005, 
van Dokkum 2005) . Our results suggest that LOSVD asym- 
metries provide fossil evidence that the final steps of the 
merger-driven galaxy growth consisted of 'dry' mergers of 
ellipticals. Previous N-body work had already noted that 
low-angular momentum E-fE mergers provide a clean way 
to reproduce many of the properties of giant, boxy ellipti- 
cals (Gonzalez-Garcia & van Albada 2005a, 2005b; Naab, 
Khochfar & Burkert 2006). Such origin for boxy ellipticals 
should be preferred over that involving bulge-less S-l-S merg- 
ers; the latter also lead to boxy remnants, but provide a poor 
match to many properties of ellipticals, including low cen- 
tral densities, too high axis ratios, and too high triaxiality 
(GGB05). 

A succession of intermediate- mass mergers, as proposed 
by Bournaud et al. (2005) for the formation of ellipticals, 
should also be investigated as a mechanism to build rem- 
nants dominated by short-axis tubes. 

Gas has been omitted from the analysis. As shown here, 
coUisionless mergers can reproduce the most salient features 
of the h^-V/a distribution without recourse to star forma- 
tion (Bekki & Shioya 1997, NBOl). Gas, however, must play 
a role in mergers of spirals, as most of the gas present in 
the precursors eventually returns back to the merger rem- 
nant (Hibbard & Mihos 1995). Given the high central con- 
centration of merger remnants, the formation of stable, or 
quasi-stable, rotationally-supported gaseous discs, is not out 
of the question. After star formation, such discs would con- 
tribute short-axis tube orbits, as NBOl point out, thus bring- 
ing models into agreement with observations. 
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